Comment on "Phase Diagram of the Random 
Energy Model with Higher-Order Ferromagnetic 
Term and Error Correcting Codes due to Sourlas" 

In a recent Letter, Dorlas and Wedagedera (DW) 0] 
have studied the random energy model (REM) with an 
additional p-spin ferromagnetic interaction, as a guide to 
the properties of a p-spin Ising model with both random 
spin glass and uniform ferromagnetic exchange, itself rel- 
evant to an error-correcting code Q. They showed that 
the non-glassy ferromagnetic phase, found for p = 2 ||] 
to lie between the paramagnetic and glassy ferromagnetic 
phases, is squeezed out to larger ferromagnetic exchange 
as p is increased and is eliminated in the limit of p — > oo. 
Here we note that (i) we have solved the corresponding 
problem of a spherical spin system with p-spin glass in- 
teractions and r-spin ferromagnetic interactions Q| and 
have shown that for all r > p > 2 the opposite situation 
applies, namely glassy ferromagnetism is suppressed and 
only non-glassy ferromagnetism remains, and (ii) a sim- 
ple mapping yields the results of DW and generalizations. 

The Hamiltonians for both the Ising and spherical 
models consist of a disordered and a ferromagnetic term: 
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where the Ji t ...i v are independent Gaussian random cou- 
plings of zero mean and variance p\J 2 /2N p ~ 1 , and (f>f — 1 
for Ising or -L </>f = 1 for spherical spins. The prop- 
erties of the system can be found from the free energy 
/sg(-^0 of the system with Jo = and a constrained 
magnetization M. They are obtained by minimizing the 
free energy 


f(M) = /sg(M) - -J M 


with respect to M, which means solving 
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Generally /sg(-^0 is first order in small M, diverges 
as \M\ — > 1, and is monotonically increasing in between. 
For r = 1, corresponding to an applied field h = Jo, 
/g G (M) = h, so the equilibrium magnetization increases 
monotonically with h and tends to unity as h — > oo. For 
r = 2, /sq(M) — JqM, so there is always a solution at 
M = 0, and a ferromagnetic solution appears continu- 
ously when J > /g G (0). For r > 2, the transition is to a 
magnetization Af m i n > 0, and M m ; n increases with r. 

The true strength of this method is in predicting 
the onset of glassiness: this depends on which parts of 
fsG(M) correspond to glassy solutions and so varies with 
model and with temperature. 


In the upper curve of the Figure we show / SG (M) for 
the REM above the glass transition temperature T s °; be- 
low T s ° the solution is glassy everywhere. At the tem- 
perature shown, the ferromagnetic transition is to a non- 
glassy phase for small enough r, while for larger r M m ; n 
is already in the glassy region. As r — > oo, J M r-1 
approaches a function which jumps from zero to Jo at 
M = 1, so A/ m in — > 1 and the transition is directly to the 
glassy ferromagnet. 



FIG. 1. Plots of dfsG/dM versus M for the random energy 
model with T/J — 0.7 (upper curve) and for the spherical 
p-spin model with p — 4 and T/J — 0.55 (lower curve). The 
solid lines represent non-glassy solutions, the dashed glassy. 

In the lower curve we show /g G (M) for the spheri- 
cal p-spin model slightly above its T s °; at some higher 
temperature the glassy region disappears; below T s ° the 
glassy region extends down to M = 0. At the tempera- 
ture shown, for small enough r, M m ; n lies in the lower 
non-glassy branch, so increasing Jo leads to a non-glassy 
ferromagnet, then a glassy, then back to a non-glassy. For 
some larger r the first non-glassy ferromagnet disappears, 
and for still larger r so does the glassy ferromagnet. A 
full calculation Q| shows that the second critical value is 
r = p. 

The discussion here has been of the static spinodal 
transition, but it can be easily extended to the thermo- 
dynamic transition by comparing the free energies @ of 
competing phases; DW concentrate on this latter case. 
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